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^ ■ Axially symmetric gravitating multi-skyrmion configurations are obtained using the 

O . harmonic map ansatz introduced in [1]. In particular, the effect of gravity on the energy 
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and baryon densities of the S'f/(3) non-gravitating multi-skyrmion configurations is studied 
in detail. 



^ . 1 Introduction 

! The Einstein- Skyrme model can be thought of as a nonlinear field theory which describes 
Q \ the interaction between a skyrmion (ie. a baryon) and a black hole. So far, mainly spheri- 
\ cally symmetric SU{N) gravitating skyrmions and black hole configurations with Skyrme 
hair have been obtained [2]-[4]. Recently, in [5, 6] SU{2) static axially symmetric reg- 
O ■ ular and black hole solutions have been derived numerically indicating that the energy 
O ! densities of such system depend on the coupling constant. In particular, as the coupling 
constant increases the energy (and baryon) density becomes either denser for the gravi- 
^1 fating skyrmions or sparser for the black hole solutions. In what follows we will construct 

the SU (3) field configurations approximating axially symmetric gravitating skyrmions and 
black holes (i.e. close to the solutions of the full equations) and investigate their properties. 
The SU{N) Einstein-Skyrme action reads: 
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+ -tr (K, Kn + -— tr {[K,, K,] [K^, K^]) 



-gd'x (1) 



where = d^UU^^ for /i = 0, 1, 2, 3; U is the SU{N) chiral field; g denotes the determi- 
nant of the metric; while k, e are coupling constants and G represents Newton's constant. 
In order for the finite-energy configurations to exist the Skyrme field has to go to a con- 
stant matrix at spatial infinity: f/ / as |a;^| oo. This effectively compactifies the 
three-dimensional Euclidean space into and hence implies that the field configurations 
of the Skyrme model can be considered as maps from into SU{N). 

The variation of the action (1) with respect to the metric g'^'^ leads to the Einstein 
equations 

R^.u - \g^.uR = SnGT^, (2) 
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with the stress-energy tensor given by 

T,.^ -y tr (X^X. - ^^^.X.X") -^^^{^^^ ~ l^^^ [^"' ^1 

(3) 

The starting point of our investigation is the introduction of the coordinates r, z, z on 
R^. In terms of the usual spherical coordinates r,9,(j) the Riemann sphere variable z is 
given by z = e*'^tan(^/2). In this system of coordinates the ansatz for the static axially 
symmetric isotropic metric reads 



m — l)z ,9 zim — I) ,9 
' dz^ + — ^ — r^dz^ + 



2{m + l) 



where the square-root of the determinant of the metric is of the form 

; — im\/l 2r^ 



dzdzj 



-g = 



/ 



[l + \z 



2^2 



(4) 



(5) 



while the metric profiles f, I, m are functions of (r, \z\'^). 

The Einstein-Skyrme system has a topological current which is covariantly conserved, 
yielding the topological charge [7] 
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and e^^"^ is the (constant) fully antisymmetric tensor. 

In the stereographically projected coordinates the action (1) for the metric (4) becomes 
equal to 
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while the baryon number takes the form 

B^ JtriKr [K„ K,]) dr dz dz. (9) 

Note that, the metric (4) leads to five nontrivial equations for the three functions /, m 
and /. In what follows, we consider static fields only (i.e. time independent Skyrme fields). 
Following the approach of [8] we consider linear superpositions of the stress-energy tensor 

Mi = T; + r,^ + r;-r°, a^2 = (r; + r^) , Ms ^ grr (t;: + t^) . (lo) 



Then, the Einstein equations (2) take the form 
87rGA^i = - + ---- + ---+ 
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where Ir denotes the partial derivative of the function I with respect to r and V the partial 
derivative with respect to the argument \z\'^ (similarly for other functions). Also, the 
expressions in (10) simplify to 
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Next we consider the static Einstein-Skyrme equations and construct field configurations 
approximating their axially symmetric solutions based on the harmonic map ansatz using 
the formulations of [1]. 



2 The Harmonic Map Ansatz 



The idea of the harmonic map ansatz (i.e. the generalisation of the rational map ansatz 
of Houghton et al [9]) involves the separation of the radial and angular dependence of the 
fields [1] as 

(13) 



U 



^2ih{r){P-I/N) 



-2ih{r)/N 



I + (e^^Kr) _ p 



3 



where P is a N x N hermitian projector which depends only on the angular variables {z, z) 
and h{r) is the corresponding profile function. Note that the matrix P can be thought of 
as a mapping from S"^ into CP^~^. Thus, P can be written as 

P{V) = (14) 

where y is a A'^-component complex vector (dependent on z and z). Following [10], we 
define a new operator P+, by its action on any vector v e (D^, as 

P,v = 3,v - (15) 

and then define further vectors P^v by induction: P^v — P+(P^~^i'). As shown in [11], 
these vectors have many interesting properties when V is holomorphic. These properties 
either follow directly from the definition of the P+ operator or are very easy to prove 
and they lead to many simplifications of the expressions for the action and the energy- 
momentum tensor. 

For (13) to be well-defined at the origin, the radial profile function h{r) has to satisfy 
/i(0) = TT while the boundary value f/ ^ J at r = oo requires that h{oo) = 0. As shown in 
[1], an attractive feature of the ansatz (13) is that it leads to a simple expression for the 
energy density which can be successively minimized with respect to the parameters of the 
projector P and then with respect to the shape of the profile function h{r). This procedure 
gives good approximations to multi-skyrmion field configurations in fiat space [1]. 

Since we are interested in axially symmetric field configurations in curved space the 
action possesses, in addition to the explicit |zp dependence due to the complex vector 
V, an implicit Iz]"^ dependence induced by the metric functions. Consequently, we cannot 
simply integrate the action over (as in flat space). Therefore, we consider a slightly 
modified ansatz by allowing the profile function h to depend on r and \z\'^ (shown in 
Figure 1). In fact, this ansatz (concerned with the Skyrme part) involves a generalisation 
of the rational map ansatz, the so-called improved harmonic map ansatz introduced in 
[12]. Indeed [12] presents better approximations (i.e. with lower energies) of pure SU{2) 
skyrmions and so are better approximants of the exact solutions obtained numerically by 
solving the full equations. 

The action (8) after the substitution of (13) with h — h{r, \z\'^) and having used the 
properties of the harmonic maps becomes 

S = J[lCr'-^-K'Ar,r' hi - k^B^, h" - (k' + ^M/i,^ + ^ -^3 ^) sin^ h 
1 sin^ h\ 

- — X — - — jdtardzdz (16) 
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Figure 1: The \z\'^ dependence of the Skyrme field profile function. [Recall, \z\ — tan ■ 
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The baryon number (9) coincides with the expression for the topological charge of the 
(jpN-i signia model (up to an overall profile dependent factor) since 
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The r integration can be performed although h — h{r, \z\'^) since siv? hhrdr — h{0)/2 
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and we have required that /i(0,|2;p) — n. 

The variation of (16) with respect to h gives its equation of motion: 
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where the dimensionless coordinate x — enr has been introduced. 

In addition, the left-hand side components (12) of the Einstein equations (11) become 
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In [9, 1] the expUcit forms of the holomorphic vector V have been obtained by minimizing 
the SU{2) and SUi^i) skyrmion energy densities, respectively. In particular, the low en- 
ergy field configurations of the SU{2) Skyrme model with baryon number from one up to 
seventeen and of the SU(?>) Skyrme model with baryon number from one up to six, which 
arc not the SU (2) embeddings, have been constructed. Prom these, the axisymmetric field 
configurations are the following: 



SU{2) B — 2 Skyrmion: The associated vector V has the simple form [9] 



V = (^^ 1)*; 



\P+V\ 
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Recently, this case has been studied in [5]. 
• SU{3) B — A Skyrmion: The harmonic vector V is of the form [1]: 
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(22) 



for a = 2.7191. The corresponding configuration has the shape of two tori on top 
of each other close to the equator of the sphere and its energy and baryon density 
are invariant under a rotation around the ;2-axis, which is equivalent to the SU{3) 
transformation: U A-'^UA for A = diag(e-2^°, 1, e'^'"). 

Next we use the vectors V given by (21) and (22) in order to obtain approximations to the 
gravitating SU {2) B = 2 and SU (3) P = 4 skyrmions, respectively. 
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3 Numerical Solutions 

Due to convenience of our numerical scheme we transform our expressions back to spherical 
coordinates where the metric is parametrized in terms of {x, 9, (f)) as 

ds-" = -fde + jdx'' + j xV + sin^ e dct>'' (23) 

while the Einstein equations (2) take the form 

fxx , I fee , 2/, 1 fl fl , cote/, , 1 fj, , 1 fek 



laMi = — — + H rrr H + --r— + 



f x^ f X f x^p p x'^ f 2 f I 2x^ f I 

2aM = + -— + —H^ _ 1 m| _^ _^ 1 /g 

2 m X m 2x'^ m 2x^ m? 2 m? 2x^ I Ax"^ P 2x^ P 

1 ^ 1 mx 1 mx Ix 1 TTie le cot 6 me cot 9 1$ 

2x I 2x m Ami m I 2x^ m x^ I 

, , 1 l^^ 1 lao \ la \ 3 It If) ^ / . X 

2q;M3 = -^ + -^^--^-|---f + — y + yCote. 24 
2 / 2x2 / 4x^ P 4 P 2x I I 

The dimensionless coupling constants a and Mj are defined as: a = 4:TtGk.'^ and Mj = 
Mi/e^K^ for i = 1,2,3. 

This system of coupled non-linear partial differential equations has to be solved numer- 
ically subject to the given boundary conditions. To map the infinite interval of the radial 
variable r onto the finite interval [0, 1] we introduce the variable x as 

Since the functions are symmetric under the transformation tt — e, it is sufficient to 
solve the equations for only < 9 < 'k/2 with the boundary conditions 

he{x, 9 = 7r/2) = /^(x, 9 = n/2) = meix, 9 = 7r/2) = leix, 9 = n/2) = . (26) 

To satisfy the regularity condition m{x, 9 = 0) = l{x, 9 = 0) we introduce the function 

, mix, 9) , , 

and impose the boundary conditions g{x,9 = 0) = 1, g{x = 0,9) = 1, g{x = 1,9) = 
1 and gg{x, 9/71/2) = 0. The numerical computations, based on the Newton-Raphson 
method, have been performed using the program FIDISOL [13] where the partial differential 
equations are discretized on a non-equidistant grid with typical grids sizes 130 x 30 and 
estimated relative numerical errors of order 10~^. 

It has been observed that, for a fixed complex vector V . the globally regular solutions 
of the Einstein-Skyrme model depend only on the dimensionless coupling constant a = 
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Figure 2: The dimensionless mass M/B in terms of the couphng constant a in the SU{2) 
and SU (3) case. 

AnGK^. In fact, the numerical simulations show that such field configurations exist only 
for a finite range of a: < a < ctmax where its maximal value depends on the specific 
form of V. For instance, in the SU(2) B — 2 case where V is given by (21) ctmax = 0.0338, 
while in the SU(3) S = 4 case where V is given by (22) with a = 2.7191 a^a^ = 0.0209. 
Assuming that our field configurations are very close to the proper solutions of the field 
equations (as shown in [12]) we believe that what we have seen here is also true for the 
exact (skyrmion) solutions. 

Hence we note that starting from a skyrmion in flat space a branch of gravitating 
skyrmions emerges for o; > 0. This branch terminates at the maximal value of the cou- 
pling constant ccmax, where it joints smoothly a second branch of solutions, which bends 
backwards from a = amax to a = 0. 



Figure 2 presents the dimensionless mass M/B = — (Gyr^/t/e) yj TQ^J\g\d'^rj /i? as 

a function of a for the SU(2) B = 2 and SU(3) B = 4 skyrmions. The mass decreases 
along the first branch with increasing a until it reaches its minimum at a = ccmax then 
it increases along the second branch with decreasing a and diverges as a tends to zero. 
Accordingly, we refer to the branches with lower (respectively higher) mass as the lower 
(respectively upper) branch. 

Figure 3 presents the a dependence of the metric functions at the origin /o = f{0,9) 
and k = /(O, 9) for the SU(2) S = 2 and the SU(3) 5 = 4 field configurations. Note 
that, /o and Iq decrease monotonically for the configurations along the lower branch with 
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Figure 3: The axial part of the initial metric profiles /(O, &) and m(0, &) as functions of 
the gravitational constant a in the SU{2) and SU{3) case. 



increasing a and along the upper branch with decreasing a; they approach finite values as 
a tends to zero on the upper branch. 

Figure 4 and 5 present the dimensionless energy density — Tq and the baryon density 
S° as functions of p = x sin ^ and z — x cos 9 for several values of the coupling constant 
a for the SU(2) S = 2 and SU(3) S = 4 skyrmions. In the SU(2) S = 2 case, both 
and possess a maximum on the p-axis which corresponds to a ring in the xy-p\ane; in 
accordance with the toroidal symmetry of the SU(2) B = 2 non-gravitating skyrmion. The 
radius of the torus decreases with increasing a on the lower branch. On the upper branch 
it continues to decrease with decreasing a and tends to zero as a tends to zero. The height 
of the maximum increases on the lower branch with increasing a and on the upper branch 
with decreasing a and diverges on the upper branch as a tends to zero. On the other hand, 
the SU(3) B = 4 gravitating skyrmion has the shape of a double torus as can be observed 
directly from Figure 5 (similar to the flat case). Note that, as a increases the radii of the 
tori and their distance from the xy-plane decrease along the lower branch and shrink to 
zero along the upper branch as a decreases. Simultaneously, the maxima of Tq and B^ 
increase and diverge on the upper branch as a tends to zero. This singular behaviour is 
seen in our choice of the dimensionless radial coordinate x = enr as first observed by Bizon 
and Chmaj [14]. If the rescaled coordinate x — xj^fa. is used instead, all functions are 
smooth and the rescaled energy M / B — ^/aM/ B remains finite as a tends to zero on the 
upper branch. 

Finally we note that, for the SU(2) B — 2 and SU(3) B — A gravitating skyrmions 
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Figure 4: The energy density e = Tq (left) and the baryon number density (right) are 
shown for the SU (2)5 = 2 solution. 
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Figure 5: Similar as Figure 3 for the SU(3) 5 = 4 solution. 
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the dependence on the couphng constant a follows the same pattern as for the spherically 
symmetric SU(2) and SU(3) skyrmions. 

4 Conclusions 

In this paper, we have used the improved harmonic map ansatz to look at field configu- 
rations involving gravitating skyrmions. In particular we have derived approximations of 
the SU{2) B — 2 and SU{3) B — A gravitating skyrmions. Our configurations are not 
solutions of the equations of motion but we believe that they are close to them; hence we 
hope that our results are close to what would have been seen for the solutions. 

We have found that, (qualitatively) the situation is very similar to the four spherical 
symmetric SU (3) skyrmions [4] since the equations have solutions only for a range of the 
gravitational coupling constant with the skyrmions being bound by the gravitational field. 
That is, for values of the coupling constant below its critical value - we have two solutions 
(of which one has much higher energy) and when the coupling constant goes beyond its 
critical value - there arc not unitary solutions for the Skyrmc fields suggesting that the 
system possesses two complex solutions. The shapes of the energy densities are very similar 
to what is seen for the non-gravitating skyrmions. Thus we suspect that the effects of the 
gravitational field are unversal in nature; the field binds the skyrmions and it alters their 
properties in a universal way making them more compact but not changing their basic 
shape. 
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